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Abstract A jet is simulated on the background of a galactic wind headed by a radiative bow shock. The wind 
shell, which is due to the radiative bow shock, is effectively destroyed by the impact of the jet cocoon, thanks 
to Rayleigh-Taylor instabilities. Associated strong HI absorption, and possibly also molecular emission, in high 
redshift radio galaxies which is observed preferentially in the smaller ones may be explained by that model, which 
is an improvement of an earlier radiative bow shock model. The model requires temperatures of ~ 10® K in the 
proto-clusters hosting these objects, and may be tested by high resolution spectroscopy of the Lya line. The 
simulations show that - before destruction - the jet cocoon fills the wind shell entirely for a considerable time 
with intact absorption system. Therefore, radio imaging of sources smaller than the critical size should reveal the 
round central bubbles, if the model is correct. 
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1. Introduction 

Radio galaxies at high redshift (z > 2) show huge (sim¬ 
ilar to the radio size), luminous 10^^ — 10"^^ erg/s) 
Lyg hal os llRottgering et al.l Il99fit Ivan Oiik et al IS 
iReuland et alJ 1200.811 Bei ng the progenitors of today’s 
brightest cluster galaxies llCarilli et al.ll200 lll. they pin¬ 
point proto-clusters of galaxies dVenemans et a,lJ 120021 
l2nn.‘ih . Many of the systems smaller than « 50 kpc show 
associated absor ption, preferentially on the blue wing of 
the emission line (Ivan O iik et al.ll997l:lBinette et al.l2000t 

iDe Breuck et al ]l2000ll . While the low column density ab¬ 

sorbers (A^hi ~ 10^^ — 10^^ cm“^) are probably mainly 
due to the Lya forest, the high column density ones 
(iVni « 10^® —10^® cm“^) are found much more frequently 
than expected from the Lya forest and hence belong to th e 
radio galaxy I van Oiik et "^119971: IWilman et al.il2003l . 
High resolution spectroscopy has revealed that some of the 
latter are made up of overla pping smaller ones. But so me 
others have been confirmed dWilman et alJl^OO.lL l2004ll . 

The frequency of this phenomenon points to common 
circumstances in the surroundings of these most massive 
objects in the young universe. A sound understanding of 
the emission line structure therefore promises insight into 
their assembly. Here, I present hydrodynamic simulations 
of a jet born inside a galactic wind shell, the latter having 
a radiative bow shock. This can be regarded as a modi¬ 


fication of an earlier model l Krausell2002ll . where the jet 
itself produced the radiative bow shock. The big absorp¬ 
tion systems might be identified with such a geometrically 
thin, dense shell, as I will discuss in the following. 

Section 13 reviews the available models from the litera¬ 
ture. Section 13 introduces the new one. The computation 
is presented in sectional and the implications and obser¬ 
vational tests are discussed in section |3 

2. Review of absorber models 

2.1. lO^^clouds in an ionisation cone 


Send offprint requests to: Martin Krause 


In this model, proposed bv ivan Oiik et al ] lll997ll . w 10^2 
clouds with a density of ric « 100 cm“®, about 10 ®Mq in 
total, are distributed throughout the region occupied by 
the radio galaxy. A central quasar, hidden from the ob¬ 
server by a central gas and dust torus, ionises the clouds 
within two cones, where the gas is also stirred up by in¬ 
teraction with the jet (observations indicate a typical line 
width of 1000 km/s). Outside the cone, the clouds are 
neither ionised by the quasar nor stirred up by the radio 
source, and hence are able to absorb the Lya emission of 
the ionised clouds in a narrower velocity range (typically 
10 — 100 km/s). The smaller sources are located in denser 
environments, which also produce more Hi absorption. 

This model cannot explain the tendency of these ab¬ 
sorption systems to be blueshifted (w 100 — 200 km/s). In 
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Figure 1. Density and temperature distribution after 80.5 Myr (0.5 Myr after jet start), 93 Myr (13 Myr after jet start) for 
the low resolution simulation. 


at least one case th e absorber has been detected in CIV 
llBinette et al ] 120001 ) ■ who inferred a lower metallicity in 
the absorber than in the emitting plasma. This would re¬ 
quire that the absorbing gas is physically distinct from the 
line emitting plasma. 

2.2. Extended low density shell 

iBinette et al ] 11200 ( 1 ) also infered a high ionisation of the 
absorbing gas. The lack of an obvious source of the ionising 
photons forced them to conclude that the shell was a low 
density quiescent absorbing screen surrounding the radio 
galaxy. They favoured a density of 10“^ — 10“^ cm“^, be¬ 
cause such an absorber would be ionised by the metagalac- 
tic background radiation. They proposed an evolutionary 
scenario: Once the radio source reaches the screen, the in¬ 
creased pressure compresses the denser parts, that from 
now on can only be seen when inside the ionisation cone, 
while the rest gets completely ionised. 

This scenario also has problems with the velocity struc¬ 
ture. Suppose there would be a quiescent gas halo at a 
temperature of a few times T = 10^ K and a density 
of n = 10 “^ cm~^ surrounding the centre of the radio 
galaxy in a distance of ~ 50 kpc. The cooling time would 


be tc = k^T/nA ss Myr, where A « 10~^^ erg cm^/ s 
is the cooling function (e.g. ISutherland fc Donital 1199,31) . 
After that time, which is much shorter than the expected 
age of the observed radio sources, a cooling flow should 
be established, and the absorbers should appear on the 
red wing, not on the blue one, as observed. It was con¬ 
cluded that the absorbing gas has to be outflowing due to 
a galactic wind ( Wilma.n et al.l2003l) . In this case, the den¬ 
sity distribution would be expected to be smooth (r:^ r~^). 
When the source expands, it should gradually reduce the 
neutral column. This has not been shown to be able to ac¬ 
commodate the observed bimodal distribution of column 
densities. Another problem arises from the temperature in 
such an outflow, which is given by 

(7 - l^)Lt 


T = 


(4/3)7rr^nfcB 
L t 


= lO'K X 


10"^^ erg/s 100 Myr n 


25 kpc 


where 7 is the adiabatic index, L is the driving power of 
the galactic wind, a few supernovae per century, r is the 
radius of the shell, and t the time for which the wind is 
active. So, one would expect the gas in these outflows to 
be too hot to contain neutral hydrogen. If they should 



























Martin Krause: Galactic Wind Shells and High Redshift Radio Galaxies 


3 





Temperature (top) and density at 93 Myr 


20 

0 

-20 



-100 -50 0 50 

Z [kpc] 

Temperature (top) and density at 97 Myr 


20 

0 

-20 


-100 -50 0 50 

Z [kpc] 


o 


o 

«• 

o 


2,74 

1.04 

-0.67 

-2.37 

-4,08 

-5,78 


o 


in 

I 


£ 

o 


«• 

o 


3.02 

1.29 

-0.43 

-2.15 

-3,87 


-5,60 


Figure 2. Density and temperature distribution after 93 Myr (13 Myr after jet start), 97 Myr (17 Myr after jet start) for the 
high resolution simulation. Only part of the grid is shown in the radial direction. Gomparison to Fig. Q shows the same overall 
evolution. Details that differ include the more evident asymmetry in the low resolution simulation, and the slower growth of 
the outer jet bubbles in the high resolution case. 


contain denser condensations, that provide the absorp¬ 
tion, the ionisation problem reappears. If the parameters 
of these winds would be tuned to the right temperature, 
the cooling flow problem reappears. The gas would be ex¬ 
pected to loose its pressure support and switch from out¬ 
flow to inflow. 


2.3. Jet with radiative bow shock 

Three consequences arise if a jet propagates into a dense 
medium. First, the bow shock slows down. Second, it may 
become radiative, an d third, very light jet physics applies 
(|k raus 3 12002112nn3|i . In this case, the bow shock is first 
spherical, following the same law of motion as the stel¬ 
lar (and galactic) wind bubble, which reads for constant 
external density p: 


.-3 


Tsh = 12 kpc 


10 "^® erg/s p 


t 


1/5 


10 Myr, 


( 1 ) 


Ush = 450km/s 


,-3 


1046 erg/s p 


/ 25 kpc 

V 1"sh 



1/3 

( 2 ) 


An external density of a few rrip cm“^ would slow the bow 
shock velocities down to the observed values. Such a slow 
bow shock would heat a surrounding X-ray atmosphere 
only marginally. The cooling time due to bremsstrahlung 
is tc = 12 MyvyjT/IhJ K/(n/cm“^). Hence, if the density 
is high enough to slow down the bow shock as required, 
the bow shock will be radiative, even if the atmosphere has 
already keV temperatures. A thin and dense shell forms, 
that may provide the required column density. Since the 
shell is subject to a number of instabilities, the shell frag¬ 
ments, and the radiation can be seen unabsorbed through 
the holes, in the larger sources. 

The natural column density of the shell, A^sh = 
10 ^^ cm“^(p/5mpcm“^)(rsh/25kpc), which would cor¬ 
respond to a surprisingly large gas mass of order 10^^ Mq, 
and is much more than observed. However, the shell 
clumps quickly, and a clumpy multi-phase medium is also 
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Figures. Golumn density of neutral hydrogen over Z-position and radial velocity (left), and integrated over space (right), 
after 80.5 Myr (top, 0.5 Myr after jet start) and 93 Myr (bottom) for the low resolution simulation. 


what is infe red from detailed simu lations of 2D radia¬ 
tive shocks I Sutherland et In fact this may be 

the only way to explain that not only atomic hydrogen 
and highly ionised carbon, bu t also molecular gas is ob - 
served at the same velocities llPe Breuck et al l lznnsaliJi . 
The effective column will be reduced. The low column re¬ 
gions might be ionised by the recently discovered inverse 
Compton (X-ray) emission of the cosmic microwave back¬ 
ground f rom the radio coco on, which may well extend into 
the UV llScharf et al.ll200,‘^ . 

Another issue is the Lya luminosity of the shell, which 
should be considerably below that of the emission line. 
Assuming one recombination per proton and a power law 
atmosphere p = this is given by: 

Ash.Lya = dTTrghUsh/lt'LyaP/TOp 


/ hz^Lya 
V "ip 


144(K-b3) 

7r^(K -|- 5)^ 


lp14 


3 X 10^3 erg/s 




25 kpc 


1046 erg/s 




rsh \ 

ro ) 
Po 


4+2k 


1/3 


5mp cm 3 


1/3 


(3) 


The observational requirements may be fulfilled, postulat¬ 
ing a high density inside of « 10 kpc, and a steep decline 


(k < —2) further out. Clearly, the atmospheres have to 
be carefully tuned . Recent estimates on jet powers easily 
reach 10^^ erg/s llChisellinil l2003t iKrause fc Camenzindl 
io^. Because high redshift radio galaxies are among the 
most powerful jet sources, this makes the model more dif¬ 
ficult. 

3. Jet within a galactic wind’s radiative bow shock 

This model involves a galactic wind due to supernova ac¬ 
tivity in the host galaxy, probably related to a starburst. 
Due to the low power of such a wind compared to a typical 
jet power, it will only have an observational effect, if the 
wind starts long before the jet activity. The wind would 
have a radiative bow shock, which is responsible for the 
aforementioned Lya absorption in the case of small jets. 
When the jets reach the shell, they destroy it so that the 
larger sources are no longer absorbed. The details of this 
process are studied in the simulation below. 

In such a model, the parameters can be well con¬ 
strained from observations. For the postshock gas to have 
a shorter cooling time than the preshock gas, the Mach 
number should be below six. This restricts the temper¬ 
ature of the atmosphere via the sound speed. A shock 
velocity of 200 km/s, implies an X-ray atmosphere with 
0.2- 1.5 X 10® K, in agreement with 
thermal bremsstrahlung in 4C 41.17 


the non-detection of 
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Column density distribution at 93 Myr 
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Figured. Golumn density of neutral hydrogen over Z-position and radial velocity (left), and integrated over space (right), 
after 93 Myr (top, 13 Myr after jet start) and 97 Myr (bottom) for the high resolution simulation. 


The age of the galactic wind bubble would be td = 
3rsh/5ush = 73Myr(rsh/25 kpc)(200 km/s/'Ush)- This is 
to exceed the cooling time due to bremsstrahlung, and 
consequently the pre-wind density is limited to n > 
0.06 cm“^y/r/10® K (25 kpc/'rsh)(t'sh/200 km/s). It could 
be even lower, if the cooling is to happen in the post¬ 
shock gas only, and more realistic cooling functions are 
applied. Such a density is similar to local clusters of 
galaxies. The resulting total absorbing column in the 
postshock gas would be two orders of magnitude below 
that of the previous model, alleviating the above de¬ 
scribed problem, but still requiring the dominant fraction 
not to be HI. The power required to drive the wind is: 
L > 5 X 10^3 gj.g/s y'T/lOe K (wsh/200 km/s)4(rsh/25 kpc), 


equivalent to some supernovae per year. Starbursts in 
otherwise normal galax ies at high redshift show winds 


with similar parameters 1 

Dawson et alJ2002l:IPettini et alJ 

l2002HTa.nken et al .112004 

). 


4. Simulation 

4.1. Setup 

The jet inside a galactic wind scenario was simu- 
lated in 2.5D w i th th e hydrodynamic code NIRVANA 
I Ziegler fc Yorkj Il997ll. To the evolution equations for 
mass, momentum and energy, the force of a dark mat¬ 


ter halo, and the zero metal, ze ro held cooling function 
from [Sutherland fc Donital lll993ll was added. The num¬ 
ber densities of electrons and ions, which determine the 
cooling rate, are computed from the total density and 
temperature, in the approximation of a pure hydrogen 
plasma in collisional ionisation equilibrium. The compu¬ 
tational domain spanned \Z x R\ = [200kpc x 100kpc], 
resolved by [2046 x 1022] cells. A control run was per¬ 
formed at double resolution, but spanned only 50 kpc 
radially, which the hrst simulation showed to be suffi¬ 
cient. The grid was initialised with an isothermal King 
atmosphere (p = 0.3 mp cm“^/(l-|-r^/lOO kpc^), = 

+ Z^) at a temperature of 10® K. In order to break the 
symmetry, random fluctuations on the percent level have 
been added to the density. The galactic wind was simu¬ 
lated by a distributed energy and mass density increase, 
(X exp(—r/3kpc), 3 x 10"^^ erg/s, and 10 Mq/jt in total. 
With these parameters, both cold gas due to a cooling 
how and a radiative bow shock by the galactic wind can 
be expected. After 80 Myr, a bipolar jet with a radius of 
1 kpc, a density of 10“® nip cin“®, a Mach number of 13, 
a velocity of 2/3 the speed of light, and a total power of 
4 X 10^® erg/s was injected in the centre of the grid. 
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Figure 5. Column over axial position (top: low resolution, dashed line: 80.5 Myr, solid line: 93 Myr; middle: high resolution, 
93 Myr; bottom: high resolution, 97 Myr). 
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Figure 6. Left: Average number density within the shell (abs(.Z)< 28 kpc) at 97 Myr for the high resolution simulation. Right: 
Maximum density, always occuring in the shell, over time. The high resolution simulation takes the higher density values. 


4.2. Results 

Density and temperature distributions at 80.5 Myr and 
93 Myr are shown in Fig. Q] (low resolution simulation) 
and at 93 Myr and 97 Myr in Fig. Ill(h igh resolution sim¬ 
ulation). 

During the first 80 Myr, the galactic wind establishes a 
radiative bow shock, which runs into a cooling flow atmo¬ 
sphere. The radiative bow shock is marginally resolved. 


i.e. the temperature rises in the shock, and then cools 
to « 10"^ K, where the cooling function was cut. Up to 
« 81 Myr, the jet stays inside of the wind shell. Then it 
quickly bores holes at the heads, and the cocoon fills up 
the wind bubble. The increased pressure inside the wind 
bubble now accelerates the shell, which is then destroyed 
due to Rayleigh-Taylor instabilities. 

In this phase, large streams of gas, raining down from 
the shell into the interior of the bubble can be observed 
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Figure 7. Covering fraction of neutral hydrogen within jZj < 30 kpc for a column in excess of 10^® cm ^ (solid), 10^® cm 
(dotted), and 10^'’ cm“^ (dashed). The left (right) panel shows the result for the low (high) resolution simulation. 
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Figures. Closeup of the shell fragments at 97 Myr for the high resolution simulation. 
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(Fig.[21and Fig.|HJ|. Being stirred up by the jet’s backflow, 
such gas may contribute to emission line gas. In the sim¬ 
ulation, even at the higher resolution, this gas soon mixes 
with the very hot and rarefied plasma from the jet cocoon, 
but in reality it may stay unspoiled at emission line tem¬ 
peratures, contributing to the emission line halo. Even so, 
the entrained gas condenses at the lower parts of the co¬ 
coon next to the jet beam (Fig.lBJ, due to the background 
gravity. The galactic wind also contributes somewhat to 
the high density in that region. Also, the gas is constantly 
stirred up due to backflowing vortices. 

The wind shell fragments into dense clumps. Hotter 
gas flows past these clumps and establishes a new, faster 
and non-radiative bow shock. A closeup of the shell frag¬ 
ments can be seen in Fig.|Hl The maximum density, which 
is reached in the shell, is shown over time in Fig. El (right). 
For the high resolution simulation, it starts to rise at 
81 Myr, when the tip of the jet hits the shell. Oscillating, 
it then keeps rising as long as the shell remains intact. In 
that phase, the jet cocoon compresses gas into the wind 
shell. At about 90 Myr, the fragmentation starts. The co¬ 
coon material can then escape from the bubble, the pres¬ 
sure drops, and the gas clouds expand again. The jet leaves 
the grid shortly after 97 Myr. Data after that time could 
be affected by that. At lower resolution, the same things 
happen, but take more time. At jet start and towards the 
end of the simulation, the maximum density is about half 
the one in the highly resolved simulation. This reflects 
the resolution increase by a factor of two. In between, the 
higher resolution simulation reaches a four times higher 
density than the low resolution one. The simulation is ob¬ 
viously not numerically converged. Nonetheless, even if 
not fully converged, the simulations serve to illustrate the 
fundamental physical behaviour that is important in un¬ 
derstanding the associated absorber phenomenon, namely 
that the jet-wind interaction causes the formation of many 
dense clouds through the operation of the R-T instability. 
Furthermore, there is good reason to apply this kind of 
resolution, which effectively cuts the R-T modes smaller 
than some critical wavelength, as explained below. 

The column density distribution for neutral hydrogen 
has been analysed at 80.5 Myr and 93 Myr (low reso¬ 
lution, Fig. EJ and 93 Myr and 97 Myr (high resolution, 
Fig.EJ. The neutral fraction has bee n calculated according 
to co llisional ionisation equilibrium ([Sutherland fc Dopital 
The 2D representations (over Z-position and R- 
velocity) show that, as long as the jet is inside the wind 
bubble, both, the cooling flow (down to -100 km/s) and 
the wind shell (200 km/s - 400 km/s) are able to form 
significant absorption systems, totally covering the radio 
galaxy. The infalling material has a very narrow veloc¬ 
ity distribution. For the wind shell, the full width at half 
maximum is w 100 km/s. The highest column density is 
found between 200 km/s and 250 km/s. The speed of the 
radiative shock is 200 km/s. At 93 Myr, the shell has been 
destroyed. At low resolution, no significant column is left 
below 500 km/s, neither of the wind shell nor of the cool¬ 
ing flow. The surviving clumps still produce high column, 


but are not coherent, neither in space nor in velocity space. 
At higher resolution, the clouds are more resilent. There 
is significant column left between 200 km/s and 400 km/s, 
also at late times. However, this comes mainly from non¬ 
central regions, and at 97 Myr the central regions are de¬ 
void of slow neutral hydrogen. 

At that time, the covering fraction is « 20%. It takes 
ss 1 — 3 Myr for a significant drop (Fig. El and Fig. □). 
At higher resolution, the covering fraction starts to drop 
later. It then drops quickly to 50%, after which it flattens, 
ending in a slightly higher number than the low resolution 
simulation. 

5. Discussion 

Using a galactic wind as the driver of the radiative shock 
reduces the lower limit on the external density to values 
observed in present day galaxy clusters. Also, the Lya self 
luminosity 0 of the shell is low, as required for absorption 
systems. The jet head pierced the wind shell after 1 Myr. 
The simulation showed that the shell is able to totally 
cover the inner 50 kpc of the jet with high column for 
7 Myr after jet start. At that time the jet is considerably 
more extended than the wind shell. This means that the 
real wind shell would be smaller than assumed here, since 
the radio jets are observed to have diameters of 50 kpc 
when the absorption disappears. 

The jet impact accelerated the absorbing shell to a 
higher velocity region, the shell thereby fragmented due 
to the Rayleigh-Taylor instability, and part of the mate¬ 
rial fell into the interior of the now radio cocoon filled 
shell. During this process the covering fraction reduced 
to w 20%. The simulations were 2D axisymmetric. In re¬ 
ality, the reduction along the third dimension should be 
similar, which would result in a final true covering fraction 
of about 20%^ w 4%. The emission lines could pass nearly 
unchanged through such a region. 

The assumption of axisymmetry also affects other 
parts of the simulation. General co nsiderations fo r jet sim¬ 
ulations are valid also in this case llNormai]lll993ll . i.e. jets 
are less stable in 3D simulations, and magnetic fields as 
well as a high Mach number may be required for stable 
propagation over the observed distances. An interesting 
extra is the interaction of the cocoon with the wind shell 
fragments. This may well cause asymmetries in the vortex 
shedding process, leading to the bending o f the jet. This is 
observed in some of the larger sou rces ('e.g. lvan Oiik et all 
Il99(it IPentericci et al.lll997l Il998ll. Therefore, further in¬ 
vestigation is necessary. 

The total neutral column is still higher than observed. 
The reason may be that a significant fraction of the 
gas cools down even further to form molecular hydrogen 
('compare IPonita fc Sutherland! l200^ p338). When ap¬ 
propriately resolved, radiative shocks have recently been 
shown to become turbu lent due to the thermal instability 
liSiitherland et al.l2003ll . Various ionisation species as well 
as molecular hydrogen might be present. Details would re¬ 
quire more involved shock models than presently available. 
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There are few examples of partially covered Lya 
regions in th e literature. TN J133 8-1942 might be a 
rare example llPe Breuck et al ] |l999t) . The higher veloc¬ 
ity clumps, present after shell destruction, may explain 
complex velocity profiles of larger sources. 

Doubling of the resolution raised the density in the 
shell by more than a factor of two. The shel l fragments are 
more resilent at higher resolution fcomnare lMellema et al.l 
l2nn2h . and therefore covering fraction and neutral column 
at low velocities reduce more slowly. The fragmentation 
of the shell starts earlier at higher resolution. This is in 
agreement with the linear growth time for the Rayleigh- 
Taylor instability, which is proportional to the wavelength, 
i.e. it starts faster if smaller wavelengths are resolved. The 
simulation is therefore not converged. Convergence would 
be achieved by processes not included in the simulation. 
The wind shell would be supported by its magnetic field, 
limiting its density rig, and the Rayleigh-Taylor instability 
would hnally be stabilised at the c ritical wavelengt h Ac by 
the transversal magnetic field B ll.Tiin et al ] |l995l) of the 
radio cocoon: 


Ac = 20 pc 


B 


10-5 G 


10 ^cm/s^\ /10cm ^ 
9 \ n-s 


■ (4) 


Here, i? is a typical shell acceleration of 500 km/s in 
10 Myr. Depending on these parameters, the shell also 
may fragment faster or slower in individual sources. For 
these reasons, the applied resolutions seem to be reason¬ 
able. 

The simulation also shows that the blueshift of the 
absorber changes with position by « 1 km/s/kpc. This 
is tiny, but might be detected in high resolution spectra. 
For sources near the critical size, high quality radio ob¬ 
servations at low frequencies may reveal the round central 
cocoon structure that is produced by the jet wind inter¬ 
action. Such observations may provide a test of the thin 
shell model. 
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